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Abstract 


We use the Miintz Legendre wavelets and operational matrix to solve a 
system of fractional integro-differential equations. In this method, the sys- 
tem of integro-differential equations shifts into the systems of the algebraic 
equation, which can be solved easily. Finally, some examples confirm- 
ing the applicability, accuracy, and efficiency of the proposed method are 
given. 
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1 Introduction 


The theory of fractional calculus and especially fractional differential equa- 
tions has recently become a popular topic, and many natural phenomena are 
modeled by it in [26, 21, 28, 11, 17, 2, 3]. The fractional integro-differential 
equations (FIDEs) are generalized integro-differential equations. In general, 
it is very difficult to obtain an exact solution for most FIDEs, so the use of 
approximate methods seems necessary. Many studies have used approximate 
methods to solve FIDE problems, such as Adomian decomposition method 
[9, 18], variational iteration method [10, 29, 22], the generalized differential 
transform method [1, 23], the homotopy perturbation method [36, 27], the 
collocation method [12, 19, 14, 34, 16], and block-pulse functions (BPFs) 
method [4, 7, 32, 33]. 

One type of the numerical methods that has been used effectively for 
solving FIDEs is wavelets methods, which can be referred to in as Legendre 
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wavelets [35], Chebyshev wavelets [37], Haar wavelets [31], Bernoulli wavelet 
[13, 30], and Miintz Legendre wavelets [5]. Consider the following nonlinear 
system of FIDEs, where s,t € [0,1],0 < p,q < 1, and D” and D® represent 
Caputo derivative: 


ee meta ee Ue ante 
D4yo(t) = for(t, yr(t), yo(t)) + fy fo2(s, u1(s), yo(s))ds 


In this article, Mintz Legendre wavelets and operational matrix have been 
used to obtain a numerical solution for relation (1). 

The organization of the paper is as follows: In Section 2, some basic 
results from the fractional calculus and the definition of Muntz Legendre 
wavelets are given. In Section 3, a numerical method based on Mintz Leg- 
endre wavelets for an approximation system of FIDEs and its convergence 
analysis is presented. In Section 4, the mentioned method has been examined 
by some examples. Finally, the conclusion is given in Section 5. 


2 Preliminaries 


In this section, some basic results from the fractional calculus and Mintz 
polynomial are given. 


2.1 Caputo derivative and Riemann—Liouville integral 


Definition 1 (see [25]). The fractional derivative of f(t) in Caputo sense for 
p,t € [0,1] is defined as 


DI) = pGqay f t-F ar 


where f(t) is an unknown function in an appropriate functional space and T 
is the gamma function. 


Definition 2 (see [25]). The Riemann-Liouville fractional integral of order 
p can be defined as follows: 


PSO) = Foy fl t-sPFle)as 


Remark 1 (see [25]). The following relationships exist between Caputo 
derivative and Riemann—Liouville integral: 
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DPI? f(t) = f(t), 


5 FPO) pn 


I? D? f(t) = f(t) — ; t>0, n-l<p<n. (2) 


! 
m=0 ue 


2.2 Muntz polynomial and Mintz Legendre wavelets 


Theorem 1 (see [24]). The sequence {t**}%o,with 0 < Ao < Al < +++ 3 00 
is fundamental in L2[0, 1] if and only if 


Co 
; eg = OW. 
k=1 


The classical Mintz polynomial is represented as 


N 
; azt*, 
k=0 


where az € R. 

The Muntz Legendre polynomial is an orthogonalized Muntz polynomial re- 
spect to the Lebesgue measure in [0, 1]. Assume that A, = {Ao, A1, A2,---,; An} 
and that 


Re(Ax) > =5 (KENo), An #Ay (RS). 


Then we can represent Mintz Legendre polynomial as follows (sce [8, 20]): 


n—-1 
: TT Ax +3 +1) 
L£y(t) = Sock, MM, Gey (n € No). 
k=0 I] (x —A;) 
J=0,j#k 


In this paper, we set A, = ky and y = 1. 
Mintz Legendre wavelets are defined on [0, 1] as (see [5]): 


7 n—1 n 
2k-1(1 + Wmy) + Lmn(2*-1t — n— 1,7), great =! < gR-7? 
0 


otherwise, 


Pnm (t) = 


where n = 1,2,...,281, m =0,1,2,...,M—1(k,M EN), and Ln(t,7) is 
the well-known Mintz Legendre of order m, that in which A; = ky. 
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3 Approximation of function by using Muntz Legendre 
wavelets 


Any function y belonging to L2[0,1] can be expanded by Mintz Legendre 
wavelets as follows (see [5]): 


[oe) Co 
y(t) = S- S- OnmPnm (Et). (3) 
n=1m=0 
Assume that 
{y1,0(t), ---;P1,—1(t), Y2,0(t), ---,92,a—1(8),---, 
(Por-1,9(t),---, Por-1, m1 (t)} 
is a set of Mintz Legendre wavelets and that 
X = span{io(t),..-, ~1,m—1(t), Y2,0(t),.-- ,%2,a—1(4),---, 
(Por-1 g(t), ---, Por-1 w_i(t)}- 


Since X is a finite-dimensional function space of L2[0,1], so y has the best 
unique approximation ym: € X such that ||y — ym/|| < ||y— || for all x € X. 
Moreover, there are unique coefficients Gnm such that 


y(t) = Ym = S> S> AnmPnm(t) = Anm®nm(t), (4) 


where 

T 
Anm = [aro, ++ +541,M—1;42,0,---,42,M—1,-++,G2k-119,--- saye-1,0—1| ) 
®nm(t) = [e10(0); 2xi9 Pie (1), Pa,0(8),24 <5 Gama b) <<: 


T 
Por-19(t),---, por-1,m—r(0)| . 


Equation (4) can be rewritten as 


where 


T 
Am’ = [eas snes anrtay ses QDMy + Gak—1(A As +9 Gn 5 
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T 
Oy = [er t),- spre (t), preg (t)s +s Pane (ts «++ Par-acargay(t)s +s mi (E) 

G=Onm, Yi = Ynm;, i=(n-1)M+m+41, m’ =2*-1 mM. 
Suppose that t; = an i = 1,2,3,...,m’ are collocation points. We define 
the Mimtz Legendre wavelets W as follows: 


1 3 5 2m’ — 1 


Vers [a )). 


2m’ 
Theorem 2 (see [6]). Suppose that y(t) € C“[0, 1]. that ym/(t) € X is the 


approximation of y(t) via wavelet, and that there exists D € N such that 
|y() (t)| < D for all ¢ € [0,1]. Then 


D 


lle, = |ly(t) — ym (OI < MioMe-p* 


Theorem 3. Suppose that in (1), fi; (4,7 = 1,2) satisfies the Lipschitz 
condition and y;_, is an approximate solution of yi(t) by the Mintz Legendre 
wavelets method. If k, M — oo, then |le;|| = ||yi — ys_,|| > 0. 


Proof. Let r;(t) (¢ = 1,2) be perturbation terms and let [p;| = m;. we can 
write y; , and y; , + e; as follows: 


m,—1 
: (k) 1 
Yi : 


| (2) filon@yin_,@))ae 


k=1 
a 7 Pi-lep 2 g bh 
+r [ (t — x) fia(s, m1, ( )ey2,( ))dsda, 5 


m,—-1 


1 : pi-l 
By fe 2 alent, (0) + e1(0),42,,(0)) + eala)de 


Pi) Jo 
1 t x pit 
Tea, [ (t— 2)" fio(s,yn_, (8) 


+ €1(8), Y2,,,(8) + €2(s))dsde. 
(7) 


By subtracting (6) from (7) and employing the Lipschitz condition, we have 


lleall <llrall + 


t 
Ly (llex|| + lleall) | (t— 2)? de 
0 


= 
(pi) 
, (8) 


t xL 
+ ——TF(|le,|| + lle eh t—2)?"dsdz. 
Tp) i({lerl] + lleall) acy (¢— 2) 
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By setting max{Z,,L2} = L, y(t t)| < Di(@t = 1,2), D = max{D},, D2} 
and using Theorem 1, we can write 


1 
lleall < |lrl] + L(lerl| 4 lle2II)(p ) 
(pi + 1) 
1 1 1 (9) 


MBMED Fp, +1) " T(pi + 2)” 


< |r; ||LD 


Since k, M — oo, then ||r;|| — 0. We can conclude that |le;|| > 0 (¢ = 
1,2). 


3.1 Operational matrix of the fractional integration 


In this section, we review the fractional-order operational matrix of integra- 
tion related to the Mintz Legendre wavelets. For this purpose, the required 
definitions and lemmas are given from [30]. 


Definition 3. For m’ = 2*-!M and i = 1,2,...,m’, the set of BPFs is 
defined as 
P a (age Ree 
b;(t) = a mo mm’ 
0 otherwise. 


We use the following properties of BPFs: 


;(2)b, (x) = to td 


0, t# J, 
: 1 _ 
[ bi(x)bj(a)dx = 4m" 
0 0, iFxj. 


Definition 4. Suppose that U = [ui,u2,..., Um]? and 
V = (v1, v2,...,Um’]?. We define 


U ® D= (wiv, U2U2,--+5 Um! Um! |" « 
Lemma 1. Let Bn = [by bay ceyBnel®, and suppose that f(t),g(t) € 
L2[0,1] can be written as f(t) = FT By (t (t Veud 9@)= GB vO’ = 
Fitter G = lan g2,0 Ga) Then 


f(a)g(a) © FT Bm (t)G? Bm (t) = (F? ® GT) Bm (t), (10) 


2 


f(a)? & (ET By (t)) = (FT PB). (11) 
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The Riemann-Liouville fractional integration of order a of BPFs can be 
presented as (see [15]): 


PP Bras (t) x FP Bm) (z), (12) 
where F” is the BPFs operational matrix with 


Lf fo-++ Sm'-1 

01 fi t+ fim'—2 
1? 1 00 1--: finr—3 
00 0 t+ fint—4 ? 


000-:-. 1 


in which fp = (k + 1)?t! — 2kPt1 + (k—1)?*1, kk =1,2,...,m’—-1. 

We now derive the operational matrix of fractional integration of Mintz 
Legendre wavelets. The integration of Mintz Legendre wavelets ®,, (a) can 
be expressed as 


t 
I®,,,/(t) = i Pn (s)ds y Pm! Pm! (t), (13) 
0 


where the m’-square matrix Pn is called Mintz Legendre wavelets opera- 
tional matrix. Also 


PD (0) & PP, Dn (l), (14) 
where P?’, is called Mintz Legendre wavelets fractional integral operational 
matrix. 


The Mintz Legendre wavelets can be expanded into m/-set BPFs as 
Bin (t) ~ Wms Bm (t). (15) 
Considering (12), we can write (14) and (15) as 


PByys (t) & TP Wyy Bry (t) = Vay I? Bry (t) & Wyn FP Bry (t), 
PP Bint (t) & LPO (t) & Wry FP Bry (t) & Wry FPUL) Oy (t). 


Finally, we conclude 
PP, Wy FOU. 
4 Examples 


In this section, some examples are solved by using the proposed method with 
different parameters. 
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Example 1. Consider 


where y1(0) = 0, y2(0) = 0. 
Exact solutions for system (16) are yi(t) = Vt and yo(t) = —Vt. Indeed 
Di y(t) & Gy @mi(t), DP ya(t) © Hy Pra (t), (17) 


where 
Gry = [91,925 93; - a 9m'|; cee = (hi, ha, hg,.. -s hm’). 


By using the initial conditions and (2), (14), (15), and (17), we have 


1 1 2 
yi(t) = 12 D2 y(t) + y1(0) © GE, P 2, Om (t) & GE, P2, Um: Bm’ (t), 
1 1 A 
yo(t) = 12 D2 yo(t) + y2(0) © HE, P?2, Om (t) © HE, P 2, Vm Bm (t). 
(18) 
Then, from (13) and (18), we obtain 
t t ‘ 1 ot 
i yi(s)ds & | GI P22, Bm (s)ds = GI, P?2, : On (s)ds 
0 0 0 (19) 
1 aL A 
at GE, P32, PL, Dyn () & GE, Pig? ini Bra’ (8), 
and similarly 
f 144 A 
| yo(s)ds © HEP.) ? Vm Bm: (t). (20) 
0 


By replacing (17)—(20) into (16), we obtain 


GE Wns Br (t) 
=TE@)[L 1.1) Ba) HGP, 
AE Win Brn (t) 


1 1 
= (GIP? Un + ALP 


m! m! 


144 


Me 2 Wal) Baw (hs 


m! 


Wm! oF Ht 


m! 


P 


Win’) Bm (t) — T(8)[1,1,...,1 
Relation (21) can be rewritten as follows: 


m’" m! m! 


1 1 
CF ge STE) (8 Yysveey Ul jag FG Pigs? Une + HEP 2 Umit, 
i 1 

AE Wm = G2,P 2, Um + HE, PZ, Vm —T(3)[1,1,... 


? 1| 1xm’" 

(22) 
By solving the nonlinear system (22), we can obtain the matrices of coeffi- 
cients G and H. 
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In Figure 1, the exact and approximate solutions using the proposed method 
(k = 6, M = 2) are plotted and the numerical values are given in Table 1. 


4 t 
—»*— Exact y, (t) 


—o— App. y(t) 


—»— Exact y,(t) 


0 0.5 1 


Figure 1: Exact and approximate solution of yi(t) and y2(t) in Example 1. 


Table 1: Exact and approximate solution and absolute error (AE) for y(t) and y2(t) 
in Example 1. 


t Ex. y(t) | Ap. y(t) | AE m(t Ex. yo(t) Ap. yo(t AE yo(t) 
2.3438e-02 | 1.5309e-01 | 1.5237e-01 | 7.2418e-04 | -1.5309e-01 | -1.5237e-01 | 7.2418e-04 
1.0156e-01 | 3.1869e-01 | 3.1861e-01 | 7.8660e-05 | -3.1869e-01 | -3.1861e-01 | 7.8660e-05 
1.7969e-01 | 4.2390e-01 | 4.2386e-01 | 3.3400e-05 | -4.2390e-01 | -4.2386e-01 | 3.3400e-05 
2.5781e-01 | 5.0775e-01 | 5.0773e-01 | 1.9431e-05 | -5.0775e-01 | -5.0773e-01 | 1.9431e-05 
3.3594e-01 | 5.7960e-01 | 5.7959e-01 | 1.3062e-05 | -5.7960e-01 | -5.7959e-01 | 1.3062e-05 
4.1406c-01 | 6.4348e-01 | 6.4347e-01 | 9.5455e-06 | -6.4348e-01 | -6.4347e-01 | 9.5455e-06 
4.9219e-01 | 7.0156e-01 | 7.0155e-01 | 7.3654e-06 | -7.0156e-01 | -7.0155e-01 | 7.3654e-06 

6 
6 
6 
6 
6 
6 


5.7031e-01 | 7.5519e-01 | 7.5518e-01 | 5.9049e-06 | -7.5519e-01 | -7.5518e-01 | 5.9049e-06 
6.4844e-01 | 8.0526c-01 | 8.0525e-01 | 4.8706e-06 | -8.0526e-01 | -8.0525e-01 | 4.8706c-06 
7.2656e-01 | 8.5239e-01 | 8.5238e-01 | 4.1065e- -8.5239e-01 | -8.5238e-01 | 4.1065e-06 
8.0469e-01 | 8.9704e-01 | 8.9704e-01 | 3.5232e- -8.9704e-01 | -8.9704e-01 | 3.5232e-06 
8.8281le-01 | 9.3958e-01 | 9.3958e-01 | 3.0660e-0 
9.7656e-01 | 9.8821e-01 | 9.8821e-01 | 2.6353e- 


-9.3958e-01 | -9.3958e-01 | 3.0660¢-06 
-9.8821e-01 | -9.8821e-01 | 2.6353¢-06 
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DP y(t) = yn (t)yo(t) + $y2?(t) + 2ya(t) — folya(s) + y2(s)]ds, 
5 


yr(t)ya(t) — yr(t) +1 — fo [yn (s) — 2y2(s)]ds, 


where y;(0) = 0, y2(0) = 0, and 0 <p,q< 1. 


(23) 


Exact solutions for the above system when p = q = 1 are y;(t) = t? and 
yo(t) =t. The exact solutions of y(t) and y2(t) for p,q € (0,1) are unknown. 
Let 


DP? y(t) © Ur ®mi(t),  Dyo(t) © Vin Omi (t); (24) 
where 
UF, = [uy, ue, ug,...,Um'], Voy = (v1, v2, U3, --- 5 Um] 


By using the initial conditions and (2), (14), (15), and (24), we have 


yi(t) = IPDP y(t) + y1(0) © UE, PE Dm (t) & UT, PP, Um Bm: (t), 
yo(t) = I7D%y2(t) + y2(0) © Vi2,P 4, Bmi(t) © VEPs Bri (t)- 


(25) 
Then, from (10), (11), (13), and (25), we obtain 
0 (Due) © UP n Brv(Q)VenPe ne B() ag 
= (U7 PU OVP Bar), 
eS 2 Qa 
y2(t) © (Vin Ph Vin Bini (t)) = (Vin Pay Unt) Bm’ (t), (27) 
t t t 
. yi(s)ds ~ | UL, PP? ®mi(s)ds =U2) PP, | ®m(s)ds 
0 0 0 (28) 
SUL Ph) SUP Va Bah, 
and in the same way 
t A 
| yo(s)ds = VP ty: Bm’ (t). (29) 
0 
By replacing (24) and (26)—(29) into (23), we obtain 
UT Um Brn (t) a gC ee Ud ® Vi PL Wm!) Brn! (t) 
+3(VEPS, Wm) Bri (t) + 2V,2,P 4, Vm Bm (t) 
sry Pint? Gm Brn (t) — Vinv Prat Um’ Brat (t), 
(30) 


Vir Vint Br (t) = 3 (Ug Phy Vm! @ ViyPs, Ut) Bm: (t) 


<UL PP? Bo Bel (t) + [A 1, Seeteng 1] ister (t) 
—UE, PUP Wn Bry (t) + 2V Pt Wy Brn (t). 


Relation (30) can be rewritten as follows: 
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UE Umi = $ (UEP Vm: @ VEVPL,Vmr) + $(VIEPS Um)” 
IVE PW yh =U LP Ca all — VP Wee 
(31) 
Vie Vine = 3 (Un Phy Vm! @ Viiv Pry Um) — Us Phyy Um’ Bm (t) 
+[1,1, 66650] ent — Veit Um + VE Pat Wm’ 


By solving the nonlinear system (31), we can obtain the matrices of co- 
efficients U and V. In Figure 2, the exact solutions for p = q = 1 and 
approximate solutions for different values of p and q are plotted. In Fig- 
ure 2, the approximate values are obtained by using the proposed method 
(k = 5, M = 3), and also the numerical values are given in Tables 2 and 3. 


1.2 1 
—x— Exact y,,(t) —*x— Exact y,(t) 
=S— p=1 o9{ O91 
—%*— p=0.85 || g=0.85 
4 |—O— p=0.7 —Og-0.7 
—— -p=0.55 0.8 | —— g=0.55 


0.7 
0.8 
0.6 
sy ry 
= 06 = 0.5 
> > 
0.4 
0.4 
0.3 
0.2 
0.2 
0.1 
of o® 
0 0.5 1 0 0.5 1 


Figure 2: Numerical solution of yi(t) and y2(t) for different values p and q in Example 
1 


Example 3. Consider 


DP y(t) = y2(t) + y3(t) — fo wr(s)ds, 
D4yo(t) = —4y3(t) — y(t) + 4 + fo yn(s)yo(s)ds, 


where yi (0) = 0, y2(0) = 1 and0<p,q<1. 


(32) 
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Table 2: Exact and approximate solution and absolute error (AE) for yi (t) in Example 
1. 


Ap.[p = 1] | AE[p = 
2.1680e-04 | 1.0829e- 
8.8949e-03 | 1.0582e-0. 
3.1460e-02 | 1.0172c-0 


Ap.|p = 0.85] | Ap.[p = 0.7] | Ap.[p = 0.55] 
9.0701e-04 3.7273e-03 1.4986e-02 
2.3452e-02 5.9316e-02 1.4077e-01 
6.8756e- 1.4291e-01 2.7310¢e-01 

6.7913e-02 | 9.5890e- 1.3193e-01 2.4150e-01 3.9923e-01 

1.1825e-01 | 8.8208e- 2.1044e-01 3.4943e-01 5.1564e-01 


t Ex.[p = | 
4 
4 
4 
5 
5 

1.8248e-01 | 7.8552e-05 3.0242e-01 4.6286e-01 6.2076e-01 
5 
5 
5 
5 
6 
5 
5 


1.0417e-02 | 1.0851e-0 
9.3750e-02 | 8.7891e- 
1.7708e-01 | 3.1359e- 
2.6042e-01 | 6.7817e- 
3.4375e-01 | 1.1816e- 
4.2708e-01 | 1.8240e- 
5.1042e-01 | 2.6053e- 
5.9375e-01 | 3.5254e- 
6.7708e-01 | 4.5844e- 
7.6042e-01 | 5.7823e-0 
8.4375e-01 | 7.1191e- 
9.2708¢c-01 | 8.5948e- 
9.8958e-01 | 9.7928¢- 


Oo 


Oo 
PRR Ree RP re RP re hd w Bl 
bw 


2.6059e-01 | 6.6793e- 4.0627e-01 5.7884e-01 7.1406e-01 
3.5259e-01 | 5.2791e- 5.2059e-01 6.9495e-01 7.9567e-01 
4.5848e-01 | 3.6398e- 6.4402e-01 8.0916e-01 8.6627e-01 
5.7825e-01 | 1.7458¢e-0 7.7524e-01 9.1976e-01 9.2690e-01 
7.1191e-01 | 4.1961e-0 9.1295e-01 1.0253e+00 9.7893e-01 
8.5945e-01 | 2.8742c- 1.0558e+00 1.1248¢+00 1.0240e+00 
9.7923e-01 | 4.9161e-05 1.1656e+00 1.1950e+00 1.0543e+00 


Table 3: Exact and approximate solution and absolute error (AE) for y2(t) in Example 
I 


t Ex[g¢g=1] | Ap.g=1] | AEl¢g=1 
1.0417e-02 | 1.0417e-02 | 1.0416e- 7.5926e- 
9.3750e-02 | 9.3750e-02 | 9.3743e-0 7.0114e-0 
1.7708e-01 | 1.7708e-01 | 1.7707e-01 | 1.3592e- 
2.6042e-01 | 2.6042e-01 | 2.6040e-01 | 2.0462c-0 
3.4375e-01 | 3.4375e-01 | 3.4372e-01 | 2.7585e- 4.1619e-01 4.7775e- 5.0415e-01 
4.2708e-01 | 4.2708e-01 | 4.2705e-01 | 3.4929e-0 4.9638e-01 5.4348e-0 5.4588e-01 


| | Ap.[g = 0.85] | Ap.[¢= 0.7] | Ap.[q = 0.55] 

7 

6 

5 

5 

5 

5 
5.1042e-01 | 5.1042c-01 | 5.1037e-01 | 4.2466ce-05 5.7259e-01 6.0188e-0 5.8140e-01 

5 

5 

5 

5 

5 

5 


2.1272¢-02 4,2951e-02 8.5180e-02 
1.4078e-01 2.0576e-0 2.8594e-01 
2.4035e-01 3.1438e- 3.8446e-01 
3.3126e-01 4.0256e-0 4.5249e-01 


LORE NO) 


5.9375e-01 | 5.9375e-01 | 5.9370e-01 | 5.0169e-0 6.4530e-01 6.5441e-0 6.1329e-01 
6.7708¢-01 | 6.7708¢e-01 | 6.7708e-01 | 5.8012c- 7.1488¢-01 7.0223¢- 6.4344e-01 
7.6042e-01 | 7.6042e-01 | 7.6035e-01 | 6.5975e-0 7.8160e-01 7.4628¢e-0 6.7336e-01 
8.4375e-01 | 8.4375c-0 8.4368¢-0 7.4034e- 8.4566c-0 7.8743¢e-0 7.0429e-01 
9.2708e-01 | 9.2708e-01 | 9.2700e-01 | 8.2170e- 9.0727e-01 8.2650e- 7.3726¢e-01 
9.8958c-01 | 9.8958c-01 | 9.8950e- 8.8310c-0 9.5195e-0 8.5492c-0 7.6385¢e-01 


Exact solutions for the above system when p = g = 1 are y;(t) = sin(t) 
and yo(t) = cos(t). The exact solutions of yi(t) and yo(t) for p,q € (0,1) are 
unknown. 

Let 
DP y(t) & Wo Omi (t), D2 yo(t) & RT, ®m(t), (33) 


where 
T T 
Wy = [W1, We, W3,---,Wm], Roy = [11,72,73,---, 1m’. 


By using the initial conditions and (2), (14), (15), and (33), we have 


yi(t) = IPD? y(t) + y1(0) & WE PP Om (t) © WE, PP, mn Bm: (t), 
yo(t) = IID %yo(t) + y2(0) & Roy PL, Om (t) +1 RE, Pt, Vm Bm (t) +1, 
(34) 
Then, from (10), (11), (13), and (34), we obtain 
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yi (t)yo(t) & (Wo, P? Win Bm: (t)) (Ro P 2, Win Bm (t) + 1) 
= (WE, PP Win @ Ro, P41, Wm) Bm (t) + Wo PP Wm Bm (t), 

(35) 

W(t) © (WE,PP Um) Bm (t), 


ya(t) © (RE/P2, Vm) Brar(t) + 2RE, PL, Win Bnr(t) + 1, 


m! 


(36) 


t t 
- yi(s)ds © | WP? Omi (s)ds & Wey PtP Un Bmi(t), (37) 
0 0 


t t 
[ edvele)as~ [WPL Wm @ Ri Phy Ym’) Brr(3)ds 
0 0) 
t 
+f W2,P? Wm (8) Bm (s)ds 
0) 
t 
= (WIP? Um: @ REP, Um) | Bm (s)ds 
0 
t A 
Wirt Pin Um! | Bm: (s)ds (38) 
0 


t 
a (Wo PE Um: @ Bi, Pt Vm) ‘| U1 Om (s)ds 
0 


t 
+ Wo PP ms | UO (s)ds 

0 
re (WEP? Wins © REPL Vyns Wat Pin Vin Bra (t) 
+ We Prt? Um Bm (t). 
By replacing (33)—(38) into (32), we obtain 


WE, Um = (WE,PP Wm) + (REPL, Wm) 
+2RF P41, Wm — Wo Pao? Um + [1,1,---51] 


m xm’? 


RE Wy = —4(RE,P2 Wy) — RE, P1, Wr — WE PP any 
—(WEPP Um @ RE, PA Wins + Wo PP Um Ut Pm Um’ 


(39) 
By solving the nonlinear system (39), we can obtain the matrices of coef- 
ficients W and R. In Figure 3, the exact solutions for p = q = 1 and 
approximate solutions for different values of p and q are plotted. In Fig- 
ure 3, the approximate values are obtained by using the proposed method 
(k = 4, M = 6), and also the numerical values are given in Tables 4 and 5. 
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Figure 3: Numerical solution of y(t), y2(t) for different values p and q in Example 2. 


Table 4: Exact and approximate solution and absolute error (AE) for yi(t) in Example 


A. 

t Ex.{p = 1] | Ap.[p=1] | AE [p=1] | Ap.[p =0.9] | Ap.[p = 0.8] | Ap.[q = 0.7] 
1.0417e-02 | 1.0416e-02 | 1.0415e-02 | 1.6952c-06 1.6784e-02 2.6928¢e-02 4.2987e-02 
9.3750e-02 | 9.3613e-02 | 9.3597e-02 | 1.5248e-05 1.2310e-01 1.6051e-01 2.0698¢e-0 
1.7708e-01 | 1.7616e-01 | 1.7613e-01 | 2.8760¢e-05 2.1677e-01 2.6388e-01 3.1639e-0 
2.6042e-01 | 2.5748e-01 | 2.5744e-0 4.2191e-05 | 3.0362e-0 3.5318e-0 4.0297e-0 
3.4375e-01 | 3.3702e-01 | 3.3696e-01 | 5.5500e-05 3.8452e-01 4.3146e-01 4.7285e-0 
4.2708e-01 | 4.1422c-01 | 4.1415e-01 | 6.8655e-05 | 4.5957e-01 4.9979e-01 5.2851e-0 
5.1042e-01 | 4.8854e-01 | 4.8846e-01 | 8.1641e-05 5.2857e-01 5.5854e-01 5.7118¢e-0 
5.9375e-01 | 5.5947e-01 | 5.5938e-01 | 9.4461e-05 | 5.9123c-01 6.0778e-01 6.0154e-0 
6.7708e-01 | 6.2652e-01 | 6.2642e-01 | 1.0715e-04 | 6.4720e-01 6.4744e-01 6.1997e-0 
7.6042e-01 | 6.8922e-01 | 6.8910e-01 | 1.1978e-04 6.9609e-01 6.7729e-01 6.2671e-0 
8.4375e-01 | 7.4714e-01 | 7.4701e-0 1.3249e-04 | 7.3749e-0 6.9709e-0 6.2197e-0 
9.2708e-01 | 7.9987e-01 | 7.9973e-01 | 1.4544e-04 7.7098e-01 7.0653e-01 6.0602e-0 
9.8958e-01 | 8.3580e-01 | 8.3564e-01 | 1.5548e-04 | 7.9065e-01 7.0661e-01 5.8697e-0 
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Table 5: Exact and approximate solution and absolute error (AE) for y2(t) in Example 


2. 

t Ex.{g=1] | Ap.{qg=1] | AE[g=1] | Ap.[q = 0.85] | Ap.[g = 0.7] | Ap.[¢ = 0.55] 
1.0417e-02 | 9.9995e-01 | 9.9989e-01 | 5.4234e-05 9.9972e-01 9.9927e- 9.9815e-01 
9.3750e-02 | 9.956le-01 | 9.9556e-01 | 5.3356e-05 9.9152e-01 9.8413¢e-0 9.7089e-01 
1.7708¢e-0 9.8436e-01 | 9.8431e-01 | 5.1034e-05 9.7364e-01 9.5657e-0 9.3023e-01 
2.6042¢-0 9.6628e-01 | 9.6624c-01 | 4.7223¢-05 9.4758e-01 9.2039e-0 8.8251e-01 
3.4375e- 9.4150e-01 | 9.4146c-01 | 4.1891e-05 9.1427e-01 8.7751¢e- 8.3059e-01 
4.2708e- 9.1018e-01 | 9.1014e-01 | 3.5016e-05 8.7448e-01 8.2939e- 7.7641e-01 
5.1042e- 8.7254e-01 | 8.7251e-01 | 2.6593e-05 8.2892e-01 7.7726e- 7.2155e-01 
5.9375e-0 8.2885e-01 | 8.2883e-01 | 1.6622e-05 7.7829e-01 7.2227e-0 6.6736e-01 
6.7708¢e- 7.7940e-01 | 7.7940e-01 | 5.1134e-06 7.2329e-01 6.6547¢- 6.1505e-01 
7.6042c- 7.2455e-0 7.2456ce-01 | 7.9178e-06 6.6465e-01 6.0790c- 5.6571¢- 
8.4375e- 6.6467e-01 | 6.6469e-01 | 2.2457e-05 6.0308e-01 5.5059e- 5.2032e-01 
9.2708e-0 6.0017e-01 | 6.0021e-01 | 3.8493e-05 5.3936e-01 4.9452e- 4.7980e-01 
9.8958e-0 5.4904e-01 | 5.4909e-01 | 5.1503e-05 4.9063e-01 4.5390e-0 4.5307e-01 


5 Conclusion 


In this paper, we introduced the Mintz Legendre wavelets to approximate the 
solution of a system of FIDEs and examined it by some numerical examples. 
Calculated absolute errors for different values of k and M indicated that by 
increasing k and M, the absolute errors decrease. The algorithm presented 
here can be easily used for different types of FIDEs. 
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